Origin of infrared peaks in the optical conductivity of ytterbium compounds 
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We have calculated optical conductivity [o"(t<j)] spectra of ytterbium compounds (YbAls, YbAl2, 
YbCu2Si2, YbNi2Ge2, YbInCu4, YbRli2Si2, YbIr2Si2, and YbBi2) based on the direct interband 
transition derived from first-principle band calculation and compared the results with the experi- 
mentally obtained spectra. The spectral feature of a peak in the middle-infrared region (mid-IR 
peak) and a shoulder structure in the far-infrared region (far-IR shoulder) in the experimental ct{ijj) 
spectra can be described by the band calculation with a common renormalization factor. This result 
indicates that the infrared spectra of Yb compounds originate from the interband transition from 
the Yb 4/ state but that the Yb 4/ state shifts to the Fermi level with strong electron correlation. 

PACS numbers: 71.27.-|~a, 79.60.-i 



In heavy fcrmion systems as well as mixed- valence ma- 
terials, such as cerium (Ce), ytterbium (Yb), and ura- 
nium compounds, a characteristic spectral feature, in 
which free carriers absorption (and consequently Drude 
weight) due to quasiparticles with heavy carrier mass and 
long relaxation time, and middle-infrared peaks (mid- 
IR peaks), appear in optical conductivity [o-(a;)] spec- 
tra in the far- and middle-infrared regions, respectivelyii 
The origin of the heavy quasiparticle state is believed to 
be the hybridization effect between conduction and local 
4/ electrons, namely c-f hybridization, which can be ex- 
plained by the periodic Anderson model (PAM) . In other 
words, the quasiparticle state due to c-f hybridization 
appears on the Fermi level (-Ep)-— The c-f hybridization 
state as well 0jS 0, heavy conduction band, which can be 
clearly explained by PAM, have been directly observed 
by a resonant angle-resolved photoemission spectroscopy 
of CeCoGei.2Sio.8, a typical heavy fermion compound.— 

On the other hand, although mid-IR peaks clearly ap- 
pear in the (t{lo) spectra of many Ce and Yb compounds,— 
their origin has remained a long-debated issue. One ex- 
planation of the origin of the mid-IR peaks is the optical 
transition between the bonding and antibonding states 
of the c-f hybridization gap described by PAM.— i^i^ In 
this case, the energy of the mid-IR peak should corre- 
spond to the hybridization intensity and be proportional 
to the square root of the Kondo temperature (Tk)i^ An- 
other explanation based on band structure has been pro- 
posed; i.e., that mid-IR peaks originate from a character- 
istic band structure in which the 4/ state is located near 
iJp.J^iiiiiS In particular, both the energy position and 
spectral shape of the mid-IR peaks of GeX^ {X = Pd, 
Sn, In) compounds can be explained by first-principle 
band structure calculation with the spin-orbit splitting 
(SOS)J^ In Yb compounds, especially YbAls, however, 
the experimental cr(w) spectra can be explained by band 
calculation in which the 4/ state is slightly shifted to- 
ward E-pt^ The origin of the energy shift of the Yb 4/ 
state is not clear, even though the energy of the Ce 4/ 
state in CeX^ is fixed for comparisonii^ 

In this paper, we describe our investigation into the 



origin of the mid-IR peak in Yb compounds and also the 
origin of the low-energy shift of the Yb 4/ state if the ori- 
gin is the band structure. We calculated ct{ijj) spectra of 
heavy fermion Yb compoimds (YbAls, YbAl2, YbCu2Si2, 
YbNi2Ge2, YbInCu4, YbRh2Si2, and YbIr2Si2) and a 
Kondo insulator (YbBi2) from first-principle band cal- 
culation with SOS and compared the results with the 
experimental a{uj) spectra. It was found that all of the 
calculated cr(w) spectra reproduce the experimental spec- 
tra as the optical transition from the Yb 4/7/2 state to 
the unoccupied state mainly originating from the Yb 5d 
state when the photon energy of all of the calculated 
spectra is reduced. This result indicates that the mid-IR 
peaks in Yb compounds also originate from the charac- 
teristic band structure. The energy shift can be regarded 
as a common renormalization factor due to the strong 
electron correlation of the Yb 4/ state. 

LDA band structure calculation was performed by the 
full potential linearized augmented plane wave plus local 
orbital (LAPW-f lo) method including SOS implemented 
in the Wien2k codeJ^ The lattice parameters listed in 
Table IJ reported in the hterature for YbAla^ YbAl2r^ 
YbCu2Si2^ YbNi2Ge2^ YbInCu4^ YbRh2Si2?i^ I- 
type YbIr2Si2i^ and YbBi^^ were used for the calcula- 
tion. RAirKmax (the smallest MT radius multiplied by 
the maximum k value in the expansion of plane waves in 
the basis set), which determines the accuracy of the basis 
set used, was set at 7.0. The total number of Brillouin 
zones was sampled with 40,000 /c-points. The obtained 
band structures of YbAlg, YbAlz, YbInCu4, YbRh2Si2, 
and YbBi2 are consistent with the previous band calcu- 
lations with SOS reported by other groups F^iiii^i^i The 
Tk is also listed in Table U for reference. Note that the 
LDA band calculation of YbBi2 produces a semi-metallic 
band structure, which is not consistent with the experi- 
mental semiconducting character. However, the purpose 
of this letter is to investigate the origin of mid-IR peaks 
in Yb compounds from the first-principle band calcula- 
tion, therefore we used the semi-metallic band structure 
from the LDA calculation. 

The (t{uj) spectra were derived from a function as 
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TABLE I: Space group, lattice constant and Kendo temperature of calculated Yb compounds. 



Space group Lattice constant (pm) Kondo temperature (K) 



YbAla 


No. 221, Pm - 3m 




a,b,c = 


420. 


670^ 


YbAla 


No. 227, Fd - 3m 




a,b,c = 


788. 


20002^ 


YbCuaSiz 


No. 139, M/mmm 


a, b 


= 391.2, 


c = 997.0i^ 


4024 


YbNiaGea 


No. 139, M/mmm 


a, b 


= 400.1, 


c = 973.3iS 


302^ 


YbInCu4 


No. 216, F - 43m 




a,b,c = 


712. 3ii 


4002S 


YbRh2Si2 No. 139, I4/mmm 


a, b 


= 400.7, 


c = 985.8i^ 


2521 


YblraSia 


No. 139, lA/mmm 


a, b 


= 403.2, 


c = 982.6i^ 


40i9 


YbBi2 


No. 225, Fm - 3m 




a,b,c = 


746.42a 


2402^ 



follows, in which direct intcrband transitions were as- 
sumedi^ 



Tie 



,{k) + iT 



'(fc) 



(1) 

Here, the \n'k) and states denote the unoccupied 
and occupied states, respectively, e and p are the po- 
larization of light and the monrentunr of the electron, 
respectively, /(e^g) is the Fermi-Dirac distribution func- 



tion, hWr, 



■ e is the energy difference between 



the unoccupied and occupied states, and V is the life- 
time. In the calculation, F = 1 meV was assumed. The 
experimental cr(ijj) spectra below the temperature of 10 K 
already reported in the literature for YbAlsf^Si YbAl2^ 
YbCu2Si2^ YbNizGcz^ YbInCu4^ YbRhzSia^ ^-type 
YbIr2Si2,— and YbBi^^ were used for comparison with 
the calculated spectra. Since the experimental cr(cj) spec- 
tra of single-crystalline YbCu2Si2, YbNi2Ge2, YbRh2Si2, 
and YbIr2Si2 have been measured in the a6-plane, (j[uj) 
calculation along the a-axis was performed. 

The calculated ct{u}) spectra of the Yb compounds 
compared with previously reported experimental spec- 
tra are shown in Fig. [TJ The experimental spectra have 
the peculiar features of the mid-IR peak (open triangles), 
as already pointed out, of the shoulder structure in the 
far-infrared region below 0.1 eV (far-IR shoulder, open 
circles) , and of the Drude weight of quasiparticlcs below 
0.02 eV. Since the calculation was only performed for 
the interband transition, the Drude weight of quasipar- 
ticlcs is not discussed here. The calculated aiu) spectra 
mainly originate from the interband transition from the 
Yb 4/7/2 to the unoccupied state including the Yb 5d 
state. The experimental and calculated horizontal scales 
are set at 0-0.5 eV (top axis) and 0-1.2 eV (bottom axis), 
respectively. From this comparison, the experimental 
(7(0;) spectra can be seen to be roughly explained by 
the calculation; for example, in YbAla, the experimental 
mid-IR peak and far-IR shoulder at around 0.25 eV and 
0.05 eV correspond to those at aroimd 0.6 eV and 0.1 eV 
in the calculation, respectively. In other Yb compounds, 
not only the mid-IR peak but also the far-IR shoulder in 
the experimental spectra are explained well by the calcu- 
lated a{u!) spectra. The experimental peak height is also 
reproduced by the band calculation except in the case of 
YbInCu4i^ For example, the heights of the mid-IR peaks 



of YbAla, YbCu2Si2, and YbBi2 exceed IQ^VL'^ cm''^ , 
which corresponds to the calculated height, while, those 
of YbRh2Si2 and YbIr2Si2 are about 5 x lO^l^-^cm-i, 
which can be reproduced by the calculation. The calcu- 
lated <j{uj) spectra can therefore successfully reproduce 
the experimental spectra. In other words, the infrared 
spectra of Yb compounds, especially the characteristic 
mid-IR peaks and far-IR shoulders, originate from the 
band structure. However, the horizontal scale of the 
calculated spectra is reduced by 0.42 (= 0.5/1.2). The 
meaning of the reducing ratio is discussed next. 

The energies of the mid-IR peaks and far-IR shoul- 
ders derived from the experimental cr(a;) spectra shown 
in Fig. [1] are plotted as a function of the calculated pho- 
ton energy of the peaks and shoulders in Fig. [2l The 
dashed line in the figure is the reducing ratio [hwexpt = 
0.42 X hojcaic) of the photon energy in the calculated a{uj) 
spectra shown in Fig. [1] This figure indicates that al- 
most all of the mid-IR peaks and far-IR shoulders except 
for the mid-IR peak of YbAl2 are located on the dashed 
line. This means that the experimental cr(a;) spectra are 
similarly renormalized from the spectra derived from the 
band structure calculation in proportion to the photon 
energy. Since both the mid-IR peaks and far-IR shoulders 
originate from the interband transition from the Yb 4/7 /2 
state, the energy shift can be explained by the renormal- 
ization of the Yb 4/7 /2 state due to the strong electron 
correlation. The Yb 4/7/2 peak in the LDA calculation 
of YbAls is actually located at about 200 mcV from the 
Ep-^, but it has been observed at 43 meV in a photoemis- 
sion spectrumi^ In YbIr2Si2 and YbRh2Si2 , the Yb 4/7 /2 
dispersion is also pushed up to Ep from that of the LDA 
calculationi^i^^i^ Therefore, since the Yb 4/7/2 peak is 
indeed strongly renormalized, the same renormalization 
effect appears in the (j{uj) spectra. 

The renormalization effect of the Yb 4/7 /2 peak can be 
described by the renormalized perturbation theory based 
on PAMj^ Due to this theory, the one-particle Green's 
function for the 4/ electron is written as 



Gfie) 



£ - e/ zA(e) - S(e) e _ g^. + iA(e) 



(2) 



where Sf and A(e) are the energy and width of the non- 
interacted 4/ level, respectively, and E(e) the self-energy 
of the 4/ state. 2, the wave function renormalization 
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FIG. 1: (Color online) Calculated optical conductivity [cr(t<j)] 
spectra compared with experimental spectra of YbAla,— 
YbAlz,-^ YbCu2Si2,-^ YbNiaGea,-^ YbInCu4,-^ YbRh2Si2?^ 
YbIr2Si2,— and YbBi2.-^ The a{u}) spectra of YbCu2Si2, 
YbNi2Ge2, YbRh2Si2, and YbIr2Si2 were calculated along the 
a-axis, in the same way as the experimental spectra. The hor- 
izontal scales of the experimental and calculated spectra are 
located at the top and bottom, respectively. Note that the 
horizontal scale of the calculated (j{u;) spectra is reduced by 
a factor of 0.42 {— 0.5/1.2) compared with that of the ex- 
perimental spectra. The vertical scales of the calculated and 
experimental a(aj) spectra are identical except in the case of 
YbInCu4, in which the experimental spectrum is three times 
as large as the calculated spectrum. The mid-IR peaks and 
far-IR shoulders of the experimental and calculated spectra 
are marked by open triangles and open circles, respectively. 



factor, is give by 



dE(e) 



de 



E=0 



(3) 



and if = zEf and A(£) = zA{e). 

Since the spectral function of the Yb 4/7/2 peak [/9/(e)] 
is derive from Pf{e) = —■^ImGf{e), the renormahzed 
Yb 4/7/2 peak (if) shifts to the low binding energy side 
from the non- interacted 4/ state (e/) by z. If the mid- 
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FIG. 2: (Color online) Experimental energies {hcjexpt) of the 
mid-IR peaks and far-IR shoulders as a function of the corre- 
sponding calculated energies {hucaic) of Yb compounds (solid 
circles). The error bars show the width of the peak top in 
both the experimental and calculated spectra. The dashed 
line indicates hujexpt = 0.42 x fiujcaic, which is the reducing 
ratio of the photon energy in Fig. [T] 



IR peaks originate from the interband transition from 
the Yb 4/7/2 to unoccupied states, £/ and if can be re- 
garded to be Hujcaic and Hujexpt, respectively. Therefore 
Hujexpt must be equal to zhiOcaic, which is consistent with 
the linear relation in Fig. [2l Then the renormalization 
factor z becomes 0.42 from the dashed line in Fig. [21 the 
self-energy E(e) = —lAe derived from Eq. [3l and the 
effective mass (m* = mg/z) becomes 2.4mo, where toq 
is the free electron mass. The renormalization effect is 
consistent with the photoemission spectrum of YbBi2, 
because the Kondo peak originating from the Yb 4/7/2 
state shifts close to E-p from the original Yb 4/7/2 state 
derived from the band structure calculation due to the 
renormalization effect4^ The same effect is considered to 
appear in the cf^lu) spectra of the other Yb compounds 
speculated from Fig. [2l This result is consistent with the 
fact that the real part of self-energy is proportional to 
the energy near and actually within a few tenths of 
an eV.— 

So far, we have clearly shown that not only the mid- 
IR peaks but also the far-IR shoulders of Yb compounds 
originate from the band structure as well as the inter- 
band transition from the Yb 4/7/2 to the unoccupied 
state. Our previous study on CeX^ {X = Pd, Sn, In) 
also pointed out that the mid-IR peaks can be explained 
by the band structure calculation^^ Therefore, it is con- 
cluded that the origin of the infrared cr(a;) peaks of both 
Ce and Yb compounds is the band structure in which the 
local 4/ state is located very near Ep. The Yb com- 
pounds need the renormalization, however, despite the 
fact that no correction is needed in the Ce compounds. 
The reason for this is not clear, but the difference in the 
4/ occupation number may play an important role. 

Next, the reason why the experimental mid-IR peak 
of YbAl2 is located at the higher energy side than the 
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dashed line in Fig.[2]is discussed. As shown in TablelH the 
Tk of YbAl2 is much higher than those of the other com- 
pounds. Both the high Tk and the intermediate valence^^ 
indicate the weak electron correlation of the Yb 4/ state 
in YbAl2. Hence, the renormalization factor is speculated 
to be close to 1. Actually, z of YbAl2 can be evaluated 
to be about 0.73 from Fig. [21 which is larger than z in 
other compounds. Therefore z is concluded to increase 
with Tk, which is consistent with the renormalized per- 
turbation theoryi^ 

The empirical common renormalization discussed in 
this work hides the physical meaning of the correlation 
effects that are generally associated with the conven- 
tional c-f hybridization and captured by several theoreti- 
cal treatments like the slave boson or the dynamical mean 
field approach. To clarify the basis of the renormaliza- 
tion, further theoretical approaches arc needed.— "^i^ 

To summarize, we derived the (j(lS) spectra of Yb com- 
pounds (YbAlg, YbAl2, YbCu2Si2, YbNizGcz, YbInCu4, 
YbRh2Si2, YbIr2Si2, and YbBi2) in the infrared region 
from first-principle band calculation and compared them 



with the experimental spectra. The calculated aiiS) spec- 
tra are consistent with the experimental spectra after the 
common renormalization even though the calculated Yb 
compounds have different crystal structures and differ- 
ent Kondo temperatures. This result indicates that the 
infrared spectra including the mid-IR peaks and far-IR 
shoulders originate from the interband transition from 
the strongly correlated Yb 4/7/2 state to the unoccupied 
state. The mid-IR peaks and far-IR shoulders reflect the 
renormalization effect of the Yb 4/7 /2 states due to the 
strong electron correlation. Although strong electronic 
correlations cannot be treated properly by LDA calcu- 
lation, they seem to be reasonable for Yb-based heavy 
fermion systems results which fit the data when the com- 
mon renormalization is applied. 
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